secret = 365658997413863400681901208936592183
p = nextprime(2^128)

{
shamir(x) = Mod(
  secret +
  40027061535687989076861383252205348682 * x +
  10380104371593008048799446356441519384 * x^2,
  p)
}

pts = [31337, 1337, 4919, 201527, 735, 13023]

{
for(i=1,#pts,
  print (pts[i], ",", lift(
    shamir(pts[i])
    ),
    ",",
    p))
}
